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The key to understanding an extrasolar giant planet's spectrum-and hence its detectability 
and evolution-lies with its atmosphere. Now that direct observations of thermal emission from 
extrasolar giant planets are in hand, atmosphere models can be used to constrain atmospheric 
composition, thermal structure, and ultimately the formation and evolution of detected planets. 
We review the important physical processes that influence the atmospheric structure and 
evolution of extrasolar giant planets and consider what has already been learned from the first 
generation of observations and modeling. We pay particular attention to the roles of cloud 
structure, metallicity, and atmospheric chemistry in affecting detectable properties through 
Spitzer Space Telescope observations of the transiting giant planets. Our review stresses the 
uncertainties that ultimately limit our ability to interpret EGP observations. Finally we will 
conclude with a look to the future as characterization of multiple individual planets in a single 
stellar system leads to the study of comparative planetary architectures. 



1. INTRODUCTION 

Atmospheres of planets serve as gatekeepers, control- 
ling the fate of incident radiation and regulating the loss 
of thermal energy. Atmospheres are also archives, preserv- 
ing gasses that reflect the formation and the evolution of 
a planet. Thus a complete characterization of an extraso- 
lar giant planet entails understanding its thermal evolution 
through time, bulk and atmospheric composition, and atmo- 
spheric structure. To date transit spectroscopy has probed 
the chemistry of the upper atmosphere of one EGP, and 
broad band measurements of the flux emitted by two ex- 
trasolar giant planets were reported in 2005. Many more 
such observations will follow as we await the direct imaging 
and resultant characterization of many EGPs around nearby 
stars. 

This review focuses on the physics of giant planet at- 
mospheres and the models which describe them. We first 
approach these planets from a theoretical perspective, pay- 
ing particular attention to those aspects of planetary models 
that directly relate to understanding detectability, character- 
ization, and evolution. We stress the modeling uncertainties 
that will ultimately limit our ability to interpret observa- 
tions. We will review the observations of the transiting giant 
planets and explore the constraints these observations place 
on their atmospheric structure, composition, and evolution. 
Unlike purely radial velocity detections, direct imaging will 
allow characterization of the atmosphere and bulk compo- 
sition of extrasolar planets, and provide data that will shed 
light on their formation and evolution through time. We will 



explore what plausibly can be learned from the first genera- 
tion of EGP observations and discuss likely degeneracies in 
interpretation that may plague early efforts at characteriza- 
tion. 

2. OVERVIEW OF GIANT PLANET ATMOSPHERES 

The core accretion theory describing the formation of gi- 
ant planets (Wetherill and Steward, 1989; Lissauer, 1993) 
suggests that any planet more massive than about 10 Earth 
masses should have accreted a gaseous envelope from the 
surrounding planetary nebula. This leads to the expecta- 
tion that any massive planet will have a thick envelope of 
roughly nebular composition surrounding a denser core of 
rock and ice. For this review we implicitly adhere to this 
viewpoint. Because subsequent processes, such as bom- 
bardment by planetesimals, can lead to enhancements of the 
heavier elements, we don't expect the composition of the 
planetary atmosphere to precisely mirror that of the nebula 
or the parent star. Observed enhancements of carbon in so- 
lar system giant planets (Figure 1), for example, range from 
a factor of about 3 at Jupiter to about 30 times solar abun- 
dance at Uranus and Neptune. 

Departures from nebular abundance provide a window 
to the formation and evolution history of a planet. The near 
uniform enrichment of heavy elements in the atmosphere 
of Jupiter (Owen et ai, 1999) has been interpreted as ev- 
idence that planetesimals bombarded the atmosphere over 
time (e.g., Atreya et ai, 2003). Direct collapse of Jupiter 
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Fig. 1. — Measured atmospheric composition of solar system giant planets (neglecting the noble gasses) expressed as a ratio to solar 
abundance (Lodders, 2003). Jupiter and Saturn abundances are as discussed in Lodders (2004), Visscher and Fegley (2005), and Flasar 
et al. (2005). Uranus and Neptune abundances are reviewed in Fegley et al. (1991) and Gautier et al. (1995). 



from nebular gas would not result in such a pattern of en- 
richment. A major goal of future observations should be 
to determine if most EGPs are similarly enriched in atmo- 
spheric heavy elements above that in the atmosphere of their 
primary stars. 

Other outstanding questions relate to the thermal struc- 
ture, evolution, cloud and haze properties, and photochem- 
istry of extrasolar giant planets. For discussion it is useful 
to distinguish between cooler, Jupiter-like planets and those 
giant planets that orbit very close to their primary stars, the 
'hot Jupiters'. While the ultimate goals for characterizing 
both types of planets are similar, the unique atmospheres 
of the two classes raise different types of questions. For 
the hot Jupiters, most research has focused on the horizon- 
tal and vertical distribution of incident stellar radiation in 
their atmospheres and the uncertain role of photochemical 
processes in altering their equilibrium atmospheric compo- 
sition. The available data from the transiting hot Jupiters 
also challenges conventional atmospheric models as their 
emergent flux seems to be grayer than expected. 

Cooler, more Jupiter-like planets have yet to be directly 
detected. Consequently most research focuses on predict- 
ing the albedos and phase curves (variation of brightness as 
a planet orbits its star caused by the angular dependence of 
atmospheric scattering) of these objects to aid their eventual 
detection and characterization. As with the solar system gi- 
ants, most of the scattered light reflecting from extrasolar 
EGPs will emerge their cloud decks. Thus developing an 
understanding of which species will be condensed at which 
orbital distances and-critically-the vertical distribution of 
those condensates is required to facilitate their characteri- 
zation. Finally, for both types of planets, second order ef- 
fects, including photochemistry and non-equilibrium chem- 
ical abundances can play surprisingly large roles in control- 
ling the observed planetary spectra. 



2.1. Atmospheric Temperature and Evolution 

A key diagnostic of the thermal state of a giant planet at- 
mosphere is the effective temperature, T c g . The total lumi- 
nosity, L, of a planet with radius R arises from the combi- 
nation of emission of absorbed incident stellar energy irJ 7 * 
and the intrinsic internal heat flux L[ nt : 

L = 4irR 2 aT^ = (1 - A)nR 2 (TrT+) + L int . (1) 

The Bond albedo, A, measures the fraction of incident en- 
ergy scattered back to space from the atmosphere. The rera- 
diation of thermalized solar photons (the first term on the 
right-hand side) makes up about 60% of the total luminos- 
ity of Jupiter and Saturn. For a hot Jupiter, 99.99% of the 
planet's luminosity is due to reradiation of absorbed stellar 
photons (the first term on the right in the equation). As a gi- 
ant planet ages, the contribution to the total luminosity from 
cooling of the interior, L lnt , falls as this energy is radiated 
away. 

For a remotely detected planet, Eq. (1) can be used to 
constrain the planetary radius. Given independent measure- 
ments of the total emitted infrared flux and the reflected 
visible flux-and assuming the internal flux is either negli- 
gible or precisely known from models-then the planetary 
radius can be inferred with an accuracy limited only by that 
of the optical and thermal IR photometry. However for 
EGPs the internal flux will not be well constrained since 
this quantity varies with the age and mass of the planet, 
both of which will not be known perfectly well. Further- 
more, for solar system planets, models only predict well 
the internal luminosity of Jupiter. Standard cooling models 
under-predict the internal luminosity of Saturn and over- 
predict that of Uranus by a sizeable margin. Realistic er- 
rors in remotely sensed EGP radii will be dominated by the 
uncertainty in L lnt and could easily exceed 25%, particu- 
larly for objects which have high internal luminosities com- 
pared to absorbed incident radiation (e.g., planets somewhat 
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younger or more massive than Jupiter). Thus bulk compo- 
sition inferred from the measured mass and radius will be 
highly uncertain. 

To aid in the interpretation of observations of a given ob- 
ject, modelers frequently compute a one-dimensional, glob- 
ally averaged temperature profile (connecting temperature 
to pressure or depth vertically through the atmosphere). Re- 
flected and emitted spectra can be computed from such a 
profile (Section 2.2). However since the fraction of the in- 
cident stellar radiation varies over a globe, one must first 
choose what fraction, /, of the normal incidence stellar con- 
stant should strike the upper layers of a one-dimensional 
atmosphere model. Setting / = 1 results in a model atmo- 
sphere that is only correct for the planet's substellar point. 
Combined with an appropriate choice for the mean solar 
incidence angle, setting / = 1/2 gives a day side average 
while / = 1/4 gives a planet- wide average. The latter is the 
usual choice for models of solar system atmospheres, since 
radiative time constants are typically long compared to ro- 
tation periods, allowing the atmosphere to come to equilib- 
rium with the mean incident flux. Such an 'average' profile 
may be less meaningful for tidally locked planets, depend- 
ing on the atmospheric temperature, radiative time constant, 
and circulation. 

The internal energy of a giant planet (Lint). a remnant of 
its formation, is transported through the bulk of the planet's 
fluid interior by efficient convection, as first discussed by 
Hubbard (1968). Whether a giant planet is at 0.05, 0.5, 
or 5 AU from its parent star, the rate at which this inter- 
nal energy is lost is controlled by the planet's atmosphere. 
In general, the closer a planet is to its parent star, or the 
smaller its flux from the interior, the deeper the boundary 
between the atmospheric radiative zone and the convective 
deep interior will be (Figure 2). Models indicate that the 
radiative/convective boundary is at ~ 0.5 bar in Jupiter and 
can range from 10 to ~ 1 kbar in hot Jupiters (Guillot et al., 
1996; Barman et al, 2001; Sudarsky et al, 2003). Cool- 
ing and contraction is slow for planets with deeper radiative 
zones because the flux carried by the atmosphere is propor- 
tional to the atmosphere's temperature gradient (see Guillot 
and Showman, 2002). 

Connecting planetary age to total luminosity or effective 
temperature presents a number of challenges. First, evo- 
lution models depend upon average planetary atmospheric 
conditions since the rate of cooling of the interior is gov- 
erned by the mean energy loss of the entire planet. For 
hot Jupiters, 'mean' conditions may involve subtleties of 
radiative transport, dynamics, and convection. Guillot and 
Showman (2002) have shown that cooling and contraction 
are hastened for models that include temperature inhomo- 
geneities at deep levels, rather than a uniform atmosphere, 
given the same incident flux. Recently, Iro et al. (2005) 
have computed time dependent radiative models for HD 
209458b, including energy transport due to constant zonal 
winds of up to 2km/s. They find that at altitudes deeper 
than the 5 bar pressure level, as the timescale for the at- 
mosphere to come into radiative equilibrium becomes very 



long, pressure-temperature profiles around the planet be- 
come uniform with longitude and time, and match a single 
'mean' profile computed using / = 1/4. This clearly sug- 
gests that model atmosphere grids computed with / = 1/4 
are most nearly correct for use as boundary conditions for 
evolution models. Models that use / = 1/2 (such as Bur- 
rows et al, 2003; Baraffe et al, 2003; Chabrier et al, 
2004), overestimate the effect of stellar irradiation on the 
evolution of giant planets, as they assume all regions of the 
atmosphere receive the flux of the day side. 

A second difficulty with evolution models relates to the 
depth at which incident stellar energy is deposited. For 
Jupiter-like atmospheres (Figure 2) the deep interior of the 
planet is connected to the visible atmosphere by a continu- 
ous adiabat. Thus absorbed stellar energy simply adds to the 
internal energy budget of the planet and-with appropriate 
book keeping-atmosphere models appropriate for isolated, 
non-irradiated objects can be used to compute the evolu- 
tion. Evolution models computed in this limit indeed work 
well for Jupiter (Hubbard, 1977). For hot Jupiters, however, 
the deep adiabat is separated by a radiative, isothermal re- 
gion from that part of the atmosphere that is in equilibrium 
with the incident radiation (Figures 2 and 3). As discussed 
by Guillot and Showman (2002), using atmosphere mod- 
els suitable for isolated objects as the boundary conditions 
for hot Jupiter evolution calculations (e.g., Burrows et al, 
2000) significantly overestimates the temperature of the at- 
mosphere at pressures of ~ 1 bar, leading one to overesti- 
mate the effect of irradiation, and predict contraction that is 
too slow. In contrast, Bodenheimer et al. (2001, 2003) and 
Laughlin et al. (2005) have computed hot Jupiter evolu- 
tion models where contraction is likely too fast. For their 
atmospheric boundary calculation, the atmospheric pres- 
sure at which the Rosseland mean optical depth reaches 2/3 
(~ 1 mbar in their models) is assigned the planetary effec- 
tive temperature, which itself is calculated after assuming a 
Bond albedo. This method assumes a very inefficient pene- 
tration of stellar flux into the planet's atmosphere, compared 
to detailed atmosphere models. Temperatures at higher at- 
mospheric pressures are underestimated, leading to an un- 
derestimation of the effect of irradation. The ideal solution 
is to compute individualized atmosphere models to use as 
boundary conditions for many timesteps in the evolutionary 
calculation. 

2.2. Spectra of Giant Planets 

The reflected and emitted spectra of giant planet atmo- 
spheres are controlled by Rayleigh and Mie scattering from 
atmospheric gases, aerosols, and cloud particles, and by 
absorption and emission of gaseous absorbers. Scattering 
of incident light usually dominates in the blue, giving way 
to absorption by the major molecular components at wave- 
lengths greater than about 0.6 /im. The major absorbers in 
the optical are methane and, for warmer planets, water. Na 
and K are major optical absorbers in the atmospheres of 
the hot Jupiters. Generally speaking, in strong molecular 
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Fig. 2. — Model temperature-pressure (T-P) profiles for cloudless atmospheres illustrating the giant planet classification scheme 
proposed by Sudarsky et al. (2003). Vertical dashed lines identify condensation curves for iron (Fe), forsterite (Mg2Si04), water 
(H2O), and ammonia (NH3) in a solar-composition atmosphere. The base of cloud for a given condensate is predicted to occur where 
the atmosphere T-P profile crosses the species' condensation curve. The Class V planets have high iron and silicate clouds, Class III 
planets are relatively cloudless, and Class I planets have high ammonia clouds, like Jupiter. 
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bands photons are absorbed before they can scatter back to 
space. In the continua between bands photons scatter be- 
fore they are absorbed. The continuum flux from a given 
object is thus controlled by Mie scattering from its clouds 
and hazes and Rayleigh scattering from the column of clear 
gas above the clouds. Longward of about 3 to 5 /^m for the 
cooler planets and at shorter wavelengths for the warmest, 
scattering gives way to thermal emission. 

The likelihood of absorption, and hence the depth of 
a given band, depends upon the molecular opacity at a 
given wavelength and the column abundance of the prin- 
cipal gaseous absorbers (e.g., methane and water) above 
a scattering layer. The column abundance of an absorber, 
in turn, depends upon the gravity (known for solar system 
planets), the height of a cloud layer, and its mixing ratio. 
Thus the spectra, even at low resolution, of EGPs are sen- 
sitive to their atmospheric temperature, metallicity, cloud 
structure, and mass. In principle by comparing observed 
spectra to models one can infer these properties from the 
data. Experience with the giant planets of our solar system, 
however, has shown that degeneracies in cloud properties 
and molecular abundances can be difficult to disentangle 
without broad, high-resolution spectral coverage. 

As with the optical and near-infrared wavelengths, the 
thermal emission of EGPs is sculpted by molecular opac- 
ities. In regions of low absorption, planets brightly emit 
from deep, warm layers of the atmosphere. Within strongly 
absorbing regions flux arises from higher, cooler layers (un- 
less there is a stratospheric temperature inversion). Bright 
emission in the window around 5 /xm was flagged by Mar- 
ley et al. (1996) as a diagnostic of the entire class of EGPs 
(but see the caveat in Sect. 3.4 below). This opacity window 
is responsible for the well-known 'five-micron hot spots' of 
Jupiter (Westphal et al., 1974). Solar system giant planets 
also exhibit true emission features arising from a tempera- 
ture inversion above the tropopause, notably in the 7.8 fim 
methane band that plays an important role in the strato- 
spheric energy budget. Photochemically produced ethane 
and acetylene also exhibit emission in some giant planet at- 
mospheres. 

The Galileo atmosphere entry probe provided a test of 
the ability of remote observers to accurately measure the 
abundance of gases in a giant planet atmosphere. Prior to 
Galileo's arrival at Jupiter the methane abundance was esti- 
mated to lie in the range of 2.0 to 3.6 times solar at Jupiter 
and 2 to 6 times solar at Saturn. Galileo measured Jupiter's 
methane abundance to be 2.9 ± 0.5 times solar (see review 
by Young 2003) and recent observations by Cassini have 
pinned Saturn's methane abundance at 10 ± 2 times solar 
(Flasar et al., 2005; Lodders, 2004). In both cases (at least 
some) remotely sensed value were accurate (e.g., Buriez 
and de Bergh, 1981). Remotely determining the abundance 
of condensed gases, such as ammonia or water, is more 
problematic and pre-Galileo measurements were not as ac- 
curate. Fortunately, ammonia will not condense in planets 
just slightly warmer than Jupiter. In young or more mas- 
sive planets, water will be in the vapor phase as well, which 



should allow for more accurate abundance retrieval. 

3. MODEL ATMOSPHERES 

Model atmospheres predict the appearance of EGPs. 
They thus facilitate the design of optimal detection strate- 
gies and play a role in interpreting observations. A typi- 
cal model recipe begins with assumptions about the atmo- 
spheric elemental composition, the atmospheric chemistry, 
the internal heat flow of the planet, the incident stellar flux, 
and various radiative transfer assumptions (e.g., is the at- 
mosphere in local thermodynamic equilibrium?). A cloud 
model for the treatment of atmospheric condensates and 
the relevant gaseous opacities are also required ingredients. 
When combined with a suitable method for handling atmo- 
spheric radiative and convective energy transport, the mod- 
eling process yields the thermal structure of the atmosphere 
and the reflected and emitted spectrum. Of course as in any 
recipe, models of giant planet atmospheres are only as good 
as the quality of the ingredients and the assumptions. Ne- 
glected physical processes, including some that might ini- 
tially seem to be of only second order importance can in fact 
have first order effects and spoil the predictions, at least in 
certain spectral regions. In this section we summarize typi- 
cal inputs into atmosphere models and discuss their relative 
contributions to the accuracy of the final product. Figure 3 
provides a comparison of hot Jupiter profiles computed by 
three different groups. The differences between the profiles 
give an indication of the uncertainty in our understanding of 
these atmospheres just due to varying modeling techniques. 

3.1. Chemistry 

Perhaps the most elemental input to an atmosphere 
model is the assumed composition of the atmosphere. To 
date most models of EGP atmospheres have assumed so- 
lar composition. However the best estimate of "solar" has 
changed over time (see Lodders, 2003 for a review) and 
of course the composition of the primary star will vary 
between each planetary system. Notably the carbon and 
oxygen abundances of the solar atmosphere remain some- 
what uncertain (Asplund, 2005) and vary widely between 
stars (Allende-Prieto et al., 2002). 

Given a set of elemental abundances, a chemical equi- 
librium calculation provides the abundances of individual 
species at any given temperature and pressure. Several sub- 
tleties enter such a calculation, particularly the treatment 
of condensates. In an atmosphere subject to a gravitational 
field, condensates are removed by settling from the atmo- 
sphere above the condensation level. Thus 'equilibrium' 
reactions that might take place were the gas to be kept 
in a sealed container will not proceed in a realistic atmo- 
sphere. A canonical example (e.g., Fegley and Lodders, 
1994) is that under pure equilibrium, sulfur-bearing gasses 
would not appear in Jupiter's atmospheres since S reacts 
with Fe grains at low temperature to form FeS. In fact, 
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Fig. 3. — Comparison of model atmosphere profiles for hot Jupiter HD209458 computed by three groups (labeled). Each profile 
assumes somewhat different intrinsic luminosity, L int , so differences at depth are not significant. All three profiles assume global 
redistribution of incident energy, or / = 4. Differences at lower pressure presumably arise from different methods for computing 
chemical equilibria, opacities, and radiative transfer assumptions. The spread in models provides an estimate of the current uncertainty 
in modeling these objects. In addition to the internal luminosity the size of the deep isothermal layer depends upon the behaviors at high 
pressure of the opacities of the major atmospheric constituents, which are poorly known. 



when iron condenses near 1600 K in Jupiter's deep atmo- 
sphere, the grains fall out of the atmosphere allowing S to 
remain as H 2 S at low temperature. The removal of conden- 
sates from the atmosphere by sedimentation is sometimes 
termed 'rainout,' but this term can be confusing since rigor- 
ously 'rain' refers only to the sedimentation of liquid water. 
Some early brown dwarf models did not properly account 
for sedimentation, but most recent modeling efforts do in- 
clude this effect (see Marley et al., 2002 for a more com- 
plete discussion). For a recent, detailed review of the atmo- 
spheric chemistry of EGPs and brown dwarfs see Lodders 
and Fegley (2006). 

Finally photochemistry, discussed further below, can al- 
ter atmospheric composition. Trace gasses produced by the 
photolysis of methane in Jupiter's atmosphere, for exam- 
ple, are both important UV absorbers and emitters in the 
thermal infrared. As such they play important roles in the 
stratospheric energy balance. Photochemical products may 
include important absorbers or hazes that may substantially 
alter the the spectra of EGPs and cloud their interpretation. 
Yelle (2004) thoroughly discusses the photochemistry of 
the upper atmosphere of EGPs and predicts thermospheres 
heated to over 10,000 K by the extreme ultraviolet flux im- 
pinging on the top of these planets' atmospheres. The high 
temperatures drive vigorous atmospheric escape by hydro- 
gen, producing an extended cloud surrounding the planet 
that has been observed in transit by Vidal-Madjar et al. 



(2003) for HD 209458 b. Despite the high escape flux, a 
negligible fraction of the total mass of the planet escapes 
over time (Yelle, 2004). 

3.2. Opacities 

For the temperature-pressure regimes found in the atmo- 
spheres of all but the hottest extrasolar giant planets, the 
most important gaseous absorbers are H 2 0, CH 4 , NH 3 , Na, 
K. In addition the pressure-induced continuum opacity aris- 
ing from collisions of H 2 with H 2 and He is particularly 
important in the thermal infrared. Other species found in 
brown dwarf atmospheres play a role in the hottest planets 
orbiting close to their primary stars. Freedman and Lod- 
ders (2006) review the current state of the various opacity 
databases used in atmospheric modeling. For the cool at- 
mospheres most likely to be directly imaged, the opacities 
are fairly well known. The greatest shortcomings of the 
current opacities are the lack of a hot molecular line list for 
CH4 and the highly uncertain treatment of the far wings of 
collisionally-broadened lines. Neither is a major limitation 
for most EGP modeling applications. 

3.3. Clouds and Hazes 

Clouds and hazes play a crucial role in controlling gi- 
ant planet spectra. In the absence of such scattering layers, 
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red photons would penetrate to deep layers of an EGP at- 
mosphere and generally be absorbed before they could be 
scattered (Marley et al, 1999), leading to very low reflec- 
tivity in the red and near-infrared. Planets with bright high 
water clouds, for example, tend to exhibit a bright contin- 
uum from scattered starlight punctuated by a few absorp- 
tion bands. Likewise silicate and iron clouds in the atmo- 
spheres of the close-in planets play major roles in control- 
ling their spectra (Seager and Sasselov, 1998). Furthermore 
for a given cloudy planet, reflected and emitted spectra are 
sensitive to the vertical distribution, fractional global cover- 
age, size distribution, and column number density of cloud 
particles. 

Unfortunately clouds are notoriously difficult to model, 
even in Earth's atmosphere where the representation of 
clouds is a leading source of uncertainty in terrestrial global 
atmospheric circulation models. Real clouds are a product 
of upward, downward, and horizontal transport of conden- 
sible vapor and solid or liquid condensate. Their detailed 
structure depends on a number of highly local factors in- 
cluding the availability of condensation nuclei and the de- 
gree of supersaturation as well as a host of microphysical 
properties of the condensate. Approaches applying a 1- 
dimensional atmosphere model to what is intrinsically a 3- 
dimensional problem are certainly overly simplistic. Never- 
theless, given the paucity of information, simple 1-D mod- 
els currently provide the most workable approach. 

A number of cloud models have been developed for so- 
lar system studies. Perhaps the most widely used has been 
an approach focusing on microphysical time constants de- 
veloped by Rossow (1978). An important shortcoming of 
such an approach is that the time constants sensitively de- 
pend upon a variety of highly uncertain factors, particularly 
the degree of supersaturation. Ackerman andMarley (2001) 
and Marley et al. (2003) review the physics employed by 
the most popular cloud models. Ackerman and Marley 
(2001) proposed a simple 1 -dimensional cloud model that 
accounts for vertical transport of condensate and conden- 
sible gas including a variable describing the efficiency of 
particle sedimentation. This model has had success fitting 
the cloudy atmospheres of L-type brown dwarfs, but is not 
able to predict such quantities as fractional global cloudi- 
ness or account for the rapidity of the L- to T-type brown 
dwarf transition. Other modeling approaches are discussed 
by Tsuji (2005), Helling et al. (2004), and Cooper et al. 
(2003) which range from purely phenomenological efforts 
to detailed numerical microphysical models of dust nucle- 
ation. Given the strong influence of clouds on EGP spec- 
tra, Sudarsky et al. (2000, 2003) have suggested classify- 
ing EGPs on the basis of which cloud layers are present 
or absent from the visible atmosphere (see Sect. 3.7). The 
suggestion is appealing but might be difficult to apply in 
practice for transitional cases, hazey planets, or for objects 
with only limited spectral data. 

Perhaps an even more challenging problem is atmo- 
spheric photochemical hazes. All of the solar system gi- 
ant planets are strongly influenced by hazes produced by 



the ultraviolet photolysis of methane. Figure 4 compares 
the incident stellar fluxes at two transiting planets with that 
received by Jupiter. The maximum wavelengths at which 
ultraviolet photons can photolyze various molecules are 
shown. At Jupiter, solar Lyman-a is an important contri- 
bution of the far UV flux. Although the primary stars of 
the hot Jupiters may lack substantial Lyman-a flux, given 
the proximity of the planets the integrated continuum radi- 
ation capable of photolyzing major molecules is compara- 
ble to or greater than that received by Jupiter. At Jupiter 
methane photolysis is the main driver of photochemistry 
since water and H 2 S are trapped in clouds far below the 
upper atmosphere. In the atmospheres of hot Jupiters, this 
will not be the case and these molecules will be rapidly 
photolyzed, perhaps providing important sources for pho- 
tochemical haze production. Ultimately haze optical depths 
depend upon production rates, condensation temperatures, 
microphysical processes, and mixing rates in the nominally 
stable stratosphere that in turn depend upon the atmospheric 
structure and poorly understood dynamical processes. To 
date only Liang et al. (2004) have considered this issue and 
then only in the context of hot Jupiter atmospheres, which 
they found to be too warm for condensation of photochem- 
ical hydrocarbon compounds. They did not consider O- 
or S-derived photochemical products. Since optically thick 
hazes can substantially alter the idealized EGP spectra and 
phase functions computed to date, much more work on their 
production is needed. In any case, disentangling the effects 
of clouds, hazes, and uncertain atmospheric abundances in 
the atmospheres of EGPs will likely require high quality 
spectra obtained over a large spectral range. 

3.4. Dynamics and Mixing 

An important limitation to conventional 1 -dimensional 
models of mean atmospheric structure is the neglect of ver- 
tical mixing. Vertical transport plays an important role 
when the dynamical time scale is short compared to a par- 
ticular chemical equilibrium timescale, as is the case for 
CO in the atmosphere of Jupiter (Prinn and Barshay, 1977; 
Fegley and Prinn, 1988; Yung et al., 1988) and cool T- 
dwarfs (Fegley and Lodders, 1996; Griffith and Yelle, 1999; 
Saumon et al, 2003; Golimowski et al, 2004). While 
methane is most abundant in the visible atmospheres, in the 
deep atmosphere, where temperatures are higher, the abun- 
dance of CO is substantially larger. Since the C-O bond is 
very strong, the conversion time from CO to CH4 in a par- 
cel of rising gas is correspondingly long. This allows ver- 
tical mixing through the atmosphere to transport CO from 
the deep atmosphere to the visible atmosphere. CO absorbs 
strongly in the M photometric band and excess CO in T 
dwarf atmospheres depresses the flux in this window region 
by up to a magnitude below that predicted by pure chemical 
equilibrium models (Saumon et al, 2003). If this mecha- 
nism also depresses the flux of cool extrasolar giant planets 
the utility of this spectral region for planet detection may 
not be as great as predicted (e.g., Burrows et al, 2001). 
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Fig. 4. — Incident flux at the top of the atmospheres of several transiting planets compared to that received by Jupiter. Vertical lines 
denote the approximate maximum wavelengths at which various molecules can be dissociated. Incident spectra at HD 149026 b and 
TrES-1 are similar to HD 209458 b and HD 189733 b (G4.5V), respectively, and are not shown for clarity. Model stellar spectra from 
Kurucz (1993). 



Another interesting effect is that of atmospheric dynam- 
ics. At lower pressures the radiative timescales are shorter 
than at higher pressures. For a tidally locked hot Jupiter, 
this will likely mean that that the upper atmosphere quickly 
adjusts to the flux it receives from the parent star, but deeper 
layers (P > 1 bar) will adjust much more sluggishly, and 
the dynamic transport of energy will be important. This is 
only beginning to be studied in detail (Showman and Guil- 
lot, 2002; Cho et al, 2003; Burkert et al, 2005; Cooper 
and Showman, 2005). Infrared observations as a function of 
planetary phase, with the Spitzer Space Telescope and per- 
haps other platforms, will enable constraints to be placed on 
atmospheric dynamics of HD 209458b and other planets. 



3.5. Albedos and Phase Curves 

Albedos are often of interest as they allow for a sim- 
ple parameterization of the expected brightness of a planet. 
Spectra of outer solar system planets with atmospheres are 
commonly reported as geometric albedo' spectra, which is 
simply the reflectivity of a planet measured at opposition. 
Other albedo definitions include the wavelength averaged 
geometric albedo and the Bond albedo, A in Eq. (1), which 
measures the ratio of scattered to incident light. Unfortu- 
nately the extrasolar planet literature on albedos has be- 
come somewhat muddled and terms are not always care- 
fully defined. Generic "albedos" are often cited with no 



definition. Yet different albedo varieties can differ from 
each other by several tenths or more. For example the com- 
monly referenced Lambert sphere has a geometric albedo of 
2/3 while an infinitely deep Rayleigh scattering atmosphere 
would have a wavelength-independent geometric albedo of 
3/4, yet both have A = 1. The two differ in the angular de- 
pendence of their scattered radiation. For absorbing atmo- 
spheres the Bond albedo depends upon the incident spec- 
trum. Since proportionately more red photons are absorbed 
than blue photons (which tend to scatter before absorption), 
an identical planet will have a different Bond albedo under 
the light of a red star rather than a blue one even though the 
geometric albedos are identical (Marley et al, 1999). 

In any case more information is needed to fully predict or 
interpret the flux observed by a distant observer of an extra- 
solar planet. Geometry dictates that extrasolar planets are 
most detectable near quadrature and not detectable at true 
opposition since they would be hidden by their star, thus 
a general description of the phase dependence of the scat- 
tered and emitted radiation is required. Phase information 
on solar system giant planets has long been used to con- 
strain cloud particle sizes and atmospheric structure. For 
example Voyager 1, which did not visit Uranus but instead 
imaged it from afar, observed the planet at high phase an- 
gles not reachable from Earth to help constrain the scatter- 
ing phase function of its clouds and Bond albedo (Pollack 
et al, 1986). Dyudina et al. (2005) recently relied upon 
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Voyager observations to derive phase curves of Jupiter and 
Saturn. Marley et al. (1999) computed phase curves for 
model EGPs by relying upon scattering tables computed by 
Dlugach and Yanovitskij (1974). More recently Sudarsky 
et al. (2005) presented a suite of model phase curves for 
EGPs. The differences between their model calculations 
and the observed phase curve Jupiter (compare their figures 
4 and 6) demonstrates that interpretation of specific plan- 
ets will always be challenging since the specifics of particle 
size and composition, hazes, and overall atmospheric struc- 
ture will likely make each giant planet discovered unique. 

3.6. Atmosphere Models 

Although pioneered by Kuiper (1952), giant planetary 
atmosphere modeling entered the modern era with the work 
of Trafton (1967) and Hogan et al. (1969). Following 
the Voyager 1 and 2 traverses of the outer solar system 
substantially more complex models were developed to ex- 
plore the atmospheric energy budgets, the thermal structure, 
and the reflected and emitted spectra (Appleby and Hogan, 
1984; Appleby, 1986; Marley and McKay, 1999) of each 
giant planet. (Note that these are forward models that com- 
bine first principle information about planetary atmospheres 
to reproduce the observed atmospheric thermal structure. 
There is also a very rich literature of inverse models that 
aid in the interpretation of specific data sets). These au- 
thors modeled the equilibrium one-dimensional radiative- 
convective thermal profiles of these atmospheres, includ- 
ing deposition of incident radiation, by modeling the atmo- 
spheric radiative transfer given observed atmospheric abun- 
dances and cloud properties. The models generally well re- 
produced observed spectra and thermal structure. This suc- 
cess provides an important reality check that 1-D modeling 
of giant planet atmospheres, given appropriate input, sat- 
isfactorily reproduces observed properties of giant planets. 
Modeling extrasolar planets, however, will be more chal- 
lenging: only the incident radiation is known with certainty. 
Atmospheric composition, cloud properties, and thermal 
structure and perhaps mass and radius will all have to be 
inferred from comparison of models to data. 

Burrows et al. (2000) reviewed the scant early work on 
atmospheric modeling of the cooler irradiated EGPs. Most 
pre- 1995 investigations focused on studying the evolution 
of isolated objects or assumed gray atmospheres to estimate 
EGP detectability. Marley (1998) computed exploratory 
spectra of irradiated giant planets and found that the pres- 
ence or absence of water clouds is an important spectral and 
albedo marker in EGP atmospheres. Particularly in the red 
and infrared the presence or absence of scattering clouds 
can change the scattered flux by a factor of two or more, 
with cloudless planets being darker. 

Atmosphere models specifically of the hot Jupiters were 
first developed by Seager and Sasselov (1998). Subsequent 
work focusing on either specific objects or the class in gen- 
eral includes that by Goukenleuque et al. (2000), Seager et 
al. (2000, 2005), Sudarsky et al. (2000, 2003), Barman et 



al. (2001), Iro et al. (2005), Burrows et al. (2005), and 
Fortney et al. (2005, 2006). As with the cooler planets, the 
main conclusion of this body of work is that the spectra of 
the hot Jupiters depends sensitively on the vertical distribu- 
tion and properties of condensates. Models that either pos- 
tulate or predict high altitude iron and silicate cloud decks 
tend to be warmer and more Planckian in thermal emission 
than models with deeper cloud decks. Hot Jupiter models 
and observations are considered in detail in Sections 3.4 and 
4. 

The most systematic surveys of model EGP spectra in- 
clude the work of Sudarsky et al. (2000, 2003, 2005) and 
Barman et al. (2001, 2005) who have studied model plan- 
ets of a variety of masses, ages, and orbital radii (Figure 
2). Burrows (2005) reviews and recasts much of the former 
work with an eye towards detectability of EGPs. The uni- 
versal conclusion of this body of work remains that molec- 
ular absorption bands and atmospheric condensates are the 
key diagnostics of giant planet effective temperature since 
giant planets cool as they age. For those planets distant 
enough from their stars that atmospheric temperature is pri- 
marily controlled by the loss of internal energy, not inci- 
dent flux, the progression to lower atmospheric temperature 
with age results in a diagnostic sequence of spectroscopic 
changes discussed in the next section. 

Planets more massive than 5 Mj may be as warm as 2000 
K shortly after formation, with temperatures falling well be- 
low 1000 K by a few hundred million years. By a few bil- 
lion years all planet mass objects (< 13Mj up it or {Burrows 
etal, 1997)) are cooler than 500 K. The important chemical 
equilibrium and condensation boundaries are shown in Fig- 
ures 2 and 3. As the atmosphere cools chemical equilibrium 
begins to favor first CH4 over CO and then NH 3 over N 2 . 
Water is present throughout this temperature range, but the 
molecular bands become deeper with falling temperature. 

The early part of this sequence has already been well 
sampled by observed L and T type brown dwarfs, the 
coolest of which is about 700 K. 

3.7. Spectral Signatures of EGPs 

No one discussion or figure can hope to capture the 
range of temperature, metallicities, and cloud structures that 
likely define the entire suite of possible giant planets. Nev- 
ertheless Figures 5 and 6 help illustrate the important phys- 
ical processes that control EGP spectra and give an indi- 
cation of the major spectral signatures expected in atmo- 
spheres with roughly solar composition. These spectra are 
purposefully presented at moderate spectral resolution that 
will likely typify early direct detection spectra. 

Figure 2 presents a set of five cloudless temperature- 
pressure (T-P) profiles. Also shown are the condensation 
curves for iron, silicate, water, and ammonia. The conden- 
sates expected in a given atmosphere depend upon the par- 
ticular atmospheric temperature structure. Sudarsky et al. 
(2003) used this atmospheric characteristic to suggest that 
planets be categorized by which clouds form in their atmo- 



9 




Fig. 5. — Model atmosphere spectra (computed by authors JF and MM) for giant planets roughly corresponding to the atmosphere 
classes shown in Figure 2. Numbers give orbital distance a from a solar type star in AU. The top model is for a hot Jupiter (a = 0.04 
AU; T e fF = 1440 K; Class V). The atmosphere is very hot with high refractory clouds. In the second model the atmosphere is cooler 
(a — 0.1 AU; T c s = 870 K; Class IV), the clouds are deeper, and the absorption bands are correspondingly more prominent. At a = 3 
AU (T ef f = 375 K; Class III) the atmosphere is cooler and relatively cloud free. Remaining two curves illustrate atmosphere with water 
clouds (a = 0.5 AU; Class II) and ammonia clouds (a = 5 AU; Class I). 



sphere. While this proposal has some drawbacks, discussed 
below, it does nicely frame the discussion of EGP atmo- 
spheres and spectra. 

The hottest EGPs orbiting most closely to their parent 
stars are expected to exhibit iron and silicate clouds high in 
their atmosphere since the atmospheric profile crosses these 
condensation curves at low pressures. Sudarsky et al. term 
such planets Class V. As seen in T-type brown dwarfs, Na 
and K are expected to be present in gaseous form and dom- 
inate the optical spectra, with water bands appearing in the 
near-infrared (Figure 4). Thermal emission is an important 
contributor to the near-infrared flux, particularly between 
the strong water bands. Cloud scattering, however, limits 
the band depths. 

In somewhat cooler atmospheres (Class IV) the clouds 
form at higher pressures in the atmosphere which results in 
deeper absorption band depths. In addition carbon is now 
found as CH 4 rather than CO; thus methane features begin 
to appear in the near infrared. 

Somewhat cooler still, the effective temperature of a 
planet orbiting a G star at 0.5 AU would be about 375 K. 
Absorption of stellar radiation keeps the atmosphere warm 
enough that water clouds would not form, yet the iron and 
silicate clouds lie far below the visible atmosphere. Al- 
though low abundance species like Na2S could form low 
optical depth hazes, these atmospheres (Class III) will be 
relatively clear with a steep blue spectral slope and deep 



molecular bands. Like Class IV and V thermal emission is 
important beyond about 1 fim.. 

For somewhat more distant planets, water and then am- 
monia condense resulting in Class II and Class I atmo- 
spheres, respectively. Continuum flux levels are controlled 
by the bright cloud decks and the 'giant planet bands' of 
methane are apparent throughout the optical, particularly 
the strong band at 0.889 /um. An ammonia absorption fea- 
ture is detectable at 1.5 yum in Class II atmospheres, but 
disappears in the colder Class I since the ammonia has con- 
densed into clouds. 

Figure 5 illustrates how sensitive such predictions are to 
atmospheric metallicity. Recall (Figure 1) that Jupiter's at- 
mosphere is enhanced by a factor of 3 in carbon and Uranus 
and Neptune by a factor of 30. The optical and near-infrared 
methane bands are highly sensitive to the methane abun- 
dance. Continuum levels, however, vary much less since 
they are controlled by the cloud structure which is not as 
sensitive to abundance variations (the clouds are already op- 
tically thick). A cloud-free atmosphere, however, is much 
darker, again illustrating the importance of clouds. 

3.8. Atmospheres of the Hot Jupiters 

Despite their nickname, "Hot Jupiters" likely bear little 
resemblance to a hotter version of our own Jupiter. Further- 
more, given their extremely small orbital separations, these 
planets have undoubtedly experienced a very different up- 
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Fig. 6. — Sensitivity of model Jupiter-like spectra (computed by authors JF and MM) to changes in model assumptions. Labeled 
curves illustrate reflected fluxes from 1 Mj planets assuming, from top to bottom, 0. 1 times solar abundance of heavy elements, solar 
abundance, 3 times solar, 5 times solar, and 1 times solar with no cloud opacity. Note that while the continuum levels are generally 
unaffected by composition changes, the depths of the methane bands are highly sensitive. The cloud-free model is much darker in the 
red and infrared since incident photons are far more likely to be absorbed at these wavelengths than Rayleigh scattered. 



bringing than their frosty jovian cousins. 

Early exploratory studies into the nature of hot-Jupiters 
revealed that stellar heating leads to much shallower T- 
P profiles than present in isolated brown dwarfs (Seager 
and Sasselov, 1998; Goukenleuque et al., 2000) (Figure 3). 
Such reductions of the temperature gradient dramatically 
weaken the strength of otherwise prominent molecular ab- 
sorption bands (e.g., due to water). Also, temperatures are 
high enough that the dominant carbon based molecule is 
CO, unlike cooler giant planet atmospheres which have high 
concentrations of CH4. Hot-Jupiter models also indicate 
that, even though significant amounts of reflected optical 
light will be present due to Rayleigh and/or Mie scattering, 
Bond albedos may be well below 0.1 (Sudarsky et al., 2000, 
2003; Barman et al, 2001). 

Even though hot-Jupiters are hot, their atmospheres are 
still cool enough that molecules, liquids and even solids 
may form and many of the issues mentioned above are still 
relevant. As always the expected atmospheric condensates 
depend on the detailed thermal strucutre of the atmosphere, 
which still varies a great deal within the hot Jupiter class. 
Equilibrium chemistry suggests that high altitude (P < 0.1 
bar) Fe and silicate clouds may be present on the day-side of 
hot-Jupiters. In general, cloud formation tends to increase 
the amount of scattered light and smooth out many of the 
spectral features. However at these altitudes the atmosphere 
is purely radiative and likely fairly quiescent which would 
allow condensate particles to quickly rain down to deeper 



levels of the atmosphere allowing only a relatively thin haze 
to remain at high altitudes. 

Hot-Jupiter temperature profiles can also enter a high 
temperature, low pressure domain in which the molecules 
TiO and VO that have strong optical absorption bands do 
not condense deeper in the atmosphere. This leads to very 
strong heating by incident radiation and the formation of ex- 
ceptionally hot stratospheres akin to, but much hotter than, 
the stratospheres driven by near-infrared methane absorp- 
tion in the solar system giant planets (Hubeny et al, 2003; 
Fortney et al., 2006). 

4. OBSERVATIONS OF HOT JUPITERS 

Since their first, surprising detection a decade ago 
(Mayor and Queloz, 1995) the hot Jupiters have received 
substantial attention, leading to the detection of the planets 
both during transit and eclipse. The chapter by Charbon- 
neau et al. fully explores this topic. Here we focus on the 
theoretical interpretation of the direct detections. 

4.1. Transmission Spectra 

As extrasolar planets transit their parent star, a small 
fraction of the stellar flux passes tangentially through the 
planet's limb and upper atmosphere. The absorbing prop- 
erties of the planetary atmosphere (along a slant geomet- 
ric path through the planet's limb) are added to the trans- 
mitted stellar absorption spectrum. There have been many 
published synthetic transmission spectra for hot-Jupiters - 
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sometimes presented as the wavelength-dependent planet 
radius that would be observed during a transit event. Some 
of these models assume plane-parallel slab geometry (Sea- 
gerandSasselov, 2000; Hubbard et al, 2001; Brown, 2001) 
while others assume spherical geometry (Barman et al, 
2001, 2002). All of these models have adopted a single 
one-dimensional thermal profile intended to represent an 
average of the planet's limb. 

Seager and Sasselov (2000) predicted strong transmis- 
sion absorption features due to Na and K alkali lines. Hub- 
bard et al. (2001) extended the modeling of transmission 
spectra to near-IR wavelengths. Their models showed that, 
similar to the Na and K alkali lines, H 2 bands can also 
imprint strong absorption feature onto the transmitted spec- 
trum. Hubbard et al. also emphasized that while the emis- 
sion spectrum of the planet may have molecular bands di- 
minished by a reduced temperature gradient, the transmis- 
sion spectrum is unaltered by such affects. 

Modeling the transmission spectrum includes many 
potential difficulties. Transmission spectroscopy probes 
the low pressure layers of the atmosphere where non- 
equilibrium conditions are most likely to occur. Also, the 
limb (or the terminator) is the transition zone between the 
night side and the irradiated day side. Consequently, the 
stellar radiation passes through a region that could have 
a steep horizontal temperature gradient and a correspond- 
ingly steep gradient in the chemical composition along tan- 
gent path lengths (see Iro et al., 2005 and Barman et al., 
2005). Such complications would be difficult to represent 
accurately using a single one-dimensional model. Fortney 
(2005) has also pointed out that, due to the relatively long 
tangential path lengths, trace condensates (negligible to the 
emission spectrum) may have a column density significant 
enough to impact the predicted transmission spectrum. 

4.2. Thermal Emission 

Impressive new datasets appeared in 2005 that placed 
new constraints on hot-Jupiter atmospheres and stress- 
tested existing hot-Jupiter atmosphere models. Those 
planets which transit their primary stars are eclipsed by 
them half an orbital period later, during which time only 
starlight-not planetary thermal emission-is detectable. The 
resulting light curve yields the ratio of planet to stellar 
flux. Observations with Spitzer Space Telescope have con- 
strained this ratio for both HD 209458b (at 24 ^m; Deming 
et al. 2005) and for TrES-1 (at 4.5 and 8 fim; Charbonneau 
et al. 2005). For HD 209458b, the known stellar flux can by 
multiplied by the flux ratio at the same wavelength to yield 
the planetary flux, 55 ± 10 /xJy, which can be expressed 
equivalently as a brightness temperature of 1130 ± 150 K. 
The distance to, and hence flux of, the TrES-1 star is not 
known. Hence only a brightness temperature can be quoted 
with some certainty. The TrES-1 temperature at 4.5 /xm is 
1010 ± 60 K and at 8 ^m is 1230 ± 60 K. 

Despite only three thermal emission data points for two 
different planets, four model interpretations (Barman et al., 



2005; Burrows et al., 2005; Fortney et al., 2005; Seager et 
al., 2005) have already been published! Some of the pub- 
lished models have conflicting interpretations. For exam- 
ple, Burrows et al. (2005) claim that the model interpreta- 
tion of TrES-1 suggests that the planet is presently reradi- 
ating on the dayside, while the best fit models of Fortney 
et al. are those for which the incident stellar radiation is 
evenly redistributed. 

Given such conflicting conclusions one might ask if the 
current data set is adequate to say anything concrete about 
the planetary atmospheres? Below we briefly summarize 
the currently published interpretations, and then provide 
our perspective on this question. The intense interest in 
HD 209458b, however, does permit comparisons between 
groups modeling the same object with different approaches. 
The range of a subset of the published models (Figure 4) 
provides a measure of the uncertainty at the current state of 
of the art. 

Burrows et al. (2005) find that their predictions for the 
planet-to-star flux density ratios of both planets to be robust 
given the uncertainties in the planets' and primary stars' 
physical properties. They inferred the presence of CO and 
perhaps H2O, and have determined that the atmospheres 
are hot. They suggest that the difference between the theo- 
retical models and all three new measurements may be ex- 
plained by an infrared-brighter hot dayside. 

Fortney et al. (2005) find that while standard solar metal- 
licity models can fit the single datapoint for HD 209458b, 
they do not for TrES-1, as the planetary spectral slope im- 
plied by the 4.5- and 8-/xm observations is redder than ex- 
pected. Model atmospheres that include a 3 to 5 times 
metal enhancement, or energy deposition into the atmo- 
sphere from 1-10 mbar, lead to a redder spectral slope. With 
these models they find they can match the TrES-1 observa- 
tions at 4.5 /im to 1 sigma, and at 8 [im to 2 sigma. Fortney 
et al. find that the best-fit models for both planets assume 
that reradiation of absorbed stellar flux occurs over the en- 
tire planet. They also note the to-date excellent agreement 
between Spitzer ultracool dwarf infrared spectral data and 
models (Roellig et al, 2004). 

In addition to standard solar abundance / = 1/4 and 
1/2 models, Barman et al. (2005) compute the 2-D vertical 
and horizontal temperature gradient over the entire day side 
in the static no-redistribution (/ = 1) case. For TrES-1, 
they find that all three models are consistent with the 8- 
/im observations at the 2 sigma level. However, only their 
/ = 1/4 model agrees with the 4.5-/im observation. More 
importantly, agreement between the / = 1 model and the 
4.5 -/zm observation would require an unrealistic change to 
the atmospheric abundances. Consequently, some process 
must be redistributing the absorbed stellar flux at a level 
more comparable to an / = 1/2 or 1/4 scenario. This is in 
agreement with the findings of Barman et al. and Fortney et 
al. (2005) for HD209458b where an/ = 1/4 model shows 
the best agreement with the 24 ztm MIPS observations. Fig- 
ure 7 also illustrates just how red TrES-1 is compared to an 
typical brown dwarf spectrum of the same luminosity. 
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Fig. 7. — The ratio of the model flux from the day side of TrES-1 to that of its parent star, from Barman et al. (2005), assuming no 
redistribution or / = 1 (top curve) and redistribution models with / = 0.5 (middle solid curve) and / = 0.25 (bottom solid curve). 
IRAC band fluxes for each model (found by convolving with the IRAC response curves) are indicated with open squares and filled circles 
show the Spitzer data with la error bars. The 4.5 pirn IRAC value for a lOx solar, / = 0.5, model is also shown (solid square). The 
lower dotted line corresponds to an isolated brown dwarf model with T e g = 1150 K. Note the much redder slope of the planet. Major 
absorption bands are indicated. 



Finally Seager et al. (2005) also conclude that a wide 
range of models fit the observational error bars. Starting 
with this philosophy, and including a 2.2- /tm observational 
upper limit {Richardson et al, 2003) neglected by the other 
groups, they rule out some models for HD 209458b at the 
hot and cold end of the plausible temperature range. They 
show that models with C/O > 1 can fit the HD 209458b 
data, including a paucity of H 2 (Figure 8) and describe 
how the same models could fit TrES-1. They suggest that 
the models show an atmospheric circulation regime inter- 
mediate between pure in situ reradiation and very efficient 
heat redistribution. 

All modelers agree on one point: hot Jupiters are in- 
deed hot. (Because the brightness temperatures of both 
HD 209458b and TrES-1 are over 1000 K, this conclusion 
does not, in fact, require a model interpretation at all.) The 
second point that all four modelers agree upon is that the 
TrES-1 4.5- and 8 -/mi data are not fit by a basic model with 
solar abundances: the model flux is too high in the 8-/mi 
band compared to observations. These two points are prob- 
ably the only concrete inferences that can be made from 
the observations. The detailed arguments for validity of 
specific atmosphere conditions (stratosphere, clouds, pho- 
tochemistry, etc.) must await further Spitzer data at other 
wavelengths. 

4.3. Non-equilibrium Effects 

So far nearly all hot-Jupiter models have assumed local 
thermodynamic equilibrium (LTE) when solving the radia- 
tive transfer equation. Assuming LTE is tantamount to as- 
suming that all species have level populations given by the 
Saha-Boltzmann distribution and that the frequency- and 
depth-dependent source function is simply a black body. 



Consequently, potentially important effects like photo- 
ionization and non-Boltzmann-like level populations are 
ignored. Given that a large fraction of the radiation field 
in a hot-Jupiter atmosphere is a non-local phenomena, LTE 
may be a rather risky assumption. Yelle (2004) reviews 
these and other issues related to understanding the likely 
very hot upper atmospheres of the hot Jupiters. 

Shortly after the detection of Na D absorption in the at- 
mosphere of HD209458b (Charbonneau et al, 2002), Bar- 
man et al. (2002) and Fortney et al. (2003) explored the 
possibility that non-LTE effects could alter the predicted 
strength of Na absorption in hot-Jupiter atmospheres. One 
possibility is that Na is ionized to pressures less than about 
0.01 bar. Also, the relevant level populations of Na may be 
underpopulated, thereby reducing Na D line strengths. 

4.4. Horizontal Gradients 

Unlike the atmospheres of isolated brown dwarfs, which 
are heated entirely from the inside out, hot Jupiters ex- 
perience significant heating by both internal (heat leftover 
from formation) and external (heat from the star) sources 
of energy. The presence of this external source of energy 
breaks the spherical symmetry implicitly assumed when 
one-dimensional model atmospheres are compared to ob- 
servations of brown dwarfs and stars, for example. A lack 
of symmetry and dual heating sources bring a number of 
challenging issues to the forefront of the hot-Jupiter model 
atmosphere problem. Indeed the day-night asymmetry and 
the need for more sophisticated modeling was recognized 
very early on, as discussed in Section 2. 1 . 

Barman et al. (2005) have made approximate 2-D static 
radiative-convective equilibrium models for the day-side at- 
mosphere for HD209458b and TrES-1. These models esti- 



13 



2500 

2000 

■S 1500 

3 1000 
E 

500 



3500 



: 


i i i | i 


t — i — i — | — i — i — i — i — | — i — i — i — i — | — i — i — i — 

LyU — i.UJ _ 






■p? ; / ' ■ W 






co i v : 




1 1 




~. i 


1 





3 3 4 

Wavelength ( M m) 



— i — i — i — i — [ — i — i — i — i — | — i — i — i — i — | — i 1 1 — i — | 1 1 — i r- 



aooo r 



c/o = 0.5 




3 3 4 

Wavelength (/J.m) 



Fig. 8. — Thermal emission spectrum for HD209458b with C/O = 1.01 (and other elements in solar abundance). For C/O > 1 at this 
planet's temperature, the water abundance is low and CH 4 abundance is increased, compared to a solar abundance spectrum (e.g., Figure 
6). For C/O greater than solar (0.5) but less than 1, the spectrum would still show reduced H2O and increased CH4 but to a lesser extent. 
The Spitzer IRAC bandpasses are shown as dotted lines. 



mate the horizontal temperature gradients-in the absence 
of winds-to be quite steep (~ 1000 K at P = 1 bar) and 
can lead to a complex chemistry gradient over the day side. 
For example, near the terminator, CO can potentially be re- 
placed as the dominate carbon bearing molecule by CH 4 . 
Na condensation may also become important in this region 
(see also Iro et al., 2005). 

As they are likely unstable, the steep horizontal gradi- 
ents in the Barman et al. (2005) model strengthen the case 
for modeling the effects of global circulation. How signif- 
icant the impact of such circulations would be on the at- 
mospheric structure (and thus on the emergent spectrum) 
depends largely on the depth at which the stellar flux is ab- 
sorb and the radiative and advective time-scales in this re- 
gion (Seager et al, 2005). 

The strong day-side irradiation has motivated several 
groups to model the global atmospheric circulation cur- 
rents in hot-Jupiter atmospheres. Three-dimensional sim- 
ulations for HD209458b suggest the possibility of strong 
zonal winds approaching or surpassing the sound speed 
(~ 1 km s _1 winds) and a significant displacement of the 
atmospheric hot spot (Showman and Guillot, 2000; Cooper 
and Showman, 2005). Additionally, Showman and Cooper 
(2005) demonstrated the impact global winds could poten- 
tially have on the depth-dependent temperature structure. 
Their simulations predict low pressure temperature inver- 
sions not seen in static 1-D atmosphere models. Cho et 
al. (2003), assuming a characteristic wind speed of 400 m 
s _1 for HD209458b, found localized jets and vortices pro- 
ducing hot and cold regions differing by as much as 300K. 
These localized "hot spots" in the Cho et al. simulations 



also move about the poles with ~ 25-day periods. Menou 
et al. (2003) applied the results of Cho et al. (2003) to other 
short-period EGPs and concluded that these kinds of circu- 
lation patterns are likely to be common among hot Jupiters. 
Despite predicting very different atmospheric flows, all of 
these simulations agree that circulation currents will most 
likely reposition the atmospheric hot spot(s) away from the 
substellar point. Consequently, maximum and minimum IR 
fluxes would not necessarily coincide with orbital phases 
that align the substellar and antistellar points with Earth, a 
result that can be tested by infrared light curves or, more 
easily, by secondary eclipse diagnostics (Williams et al., 
2006). 

While the various hydrodynamic simulations differ sig- 
nificantly in details, there appears to be agreement that vari- 
ations on the order of ~ 300 to 500 K may be present at 
"photospheric" pressures. The studies of Showman and 
Guillot (2001) and more recently Cooper and Showman 
(2005) predict steady eastward supersonic winds producing 
an atmospheric hot spot that may be displaced 60° from the 
planet's substellar point. 

5. THE FUTURE 

Transiting planets and hot Jupiters are the most favorable 
extrasolar planets for observational study in the near future. 
They can be studied in the combined light of the planet- 
star system without direct imaging. Although ground-based 
observations will surely continue, the stable conditions of 
space make it the best place for the requisite high-precision 
observations. 
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The transiting planets orbiting the brightest stars, partic- 
ularly HD209458b (the touch stone for hot Jupiters) as well 
as the newly discoved HD 189733b, will certainly continue 
to receive great attention. Spitzer photometric and spectral 
observations will determine if there are phase variations in 
the thermal emission and will search for spectral signatures 
of the atmosphere. At visible wavelengths, the Canadian 
MOST space telescope will observe HD 209458b during 
secondary eclipse in scattered light and will reach a geo- 
metric albedo of 0.15. HST data will add to the variety of 
data on the same planet. 

Other non-transiting hot Jupiters will likely be moni- 
tored with Spitzer for phase variation with both IRAC and 
MIPS. New transiting planets around bright stars, such as 
TrES-1 and HD 149026 will also be observed with Spitzer 
and HST. SOFIA and large ground-based observatories such 
as Keck and VLT may be able to detect thermal emission 
at shorter wavelengths than Spitzer (2-4 //m) during sec- 
ondary eclipse of transiting hot Jupiters. 

The more distant future for extrasolar planet characteri- 
zation looks even more promising. In the next decade Ke- 
pler will find dozens of transiting extrasolar giant planets 
at a variety of semi-major axes out to 1 AU. Some of those 
planets' atmospheres may be detectable by JWST. Further 
to the future, in the later part of the next decade, planned 20- 
to 30-m ground-based telescopes may be able to directly 
image massive Jupiters in Jupiter-like orbits. Around the 
same time, NASA's Terrestrial Planet Finders and ESA's 
Darwin aim to directly detect and characterize nearby extra- 
solar planets ranging from giant planets down to Earth-size 
planets. 

As the era of extrasolar giant planet characterization 
moves from the hot Jupiters to the realm of true Jupiter 
analogs, the focus of characterization will change. Mass 
and particularly radius, which are most easily constrained 
by transits, will be much less well constrained for more 
distant planets. Instead spectra, which are especially sen- 
sitive to cloud structure, atmospheric composition, and at- 
mospheric chemistry, will provide the primary method for 
planet characterization. Ultimately obtaining the equivalent 
of the compositional fingerprint shown in Figure 1 for many 
planets in many planetary systems will illuminate the gi- 
ant planet formation process as a function of stellar mass, 
metallicity, and planetary system architecture. 
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